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Abstract 

In present paper, stress-strain relationship of <001> oriented 
PMN-0.32PT relaxor ferroelectric single crystals has been 
investigated. Polarization rotation model is developed to 
explain the observed behaviors of PMN-0.32PT. Stress-strain 
curve of ferroelectric single crystal BaTi03 is calculated with 
the first principle method. Results show that the polarization 
rotation model is reasonable. Based on the polarization 
rotation model, a micromechanical constitutive model of 
PMN-0.32PT is constructed. It is shown that the developed 
model can faithfully capture the typical characteristic of the 
observed behavior of <001> oriented PMN-0.32PT. 
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Introduction 

Note that PZN-PT and PMN-PT single crystals usually 
experience electric and/or mechanical loading during 
their in-service life. It has been shown that externally 
applied loading has significant effect on the properties 
of these crystals[l]. It has been shown loading in the 
form of electric field[l-3], and stress[4], can lead to 
polarization switching and phase transition, which 
changes dramatically alter their electromechanical 
properties. A mature level understanding of their 
responses to electrical, mechanical and temperature 
loading condition is thus essential to fulfill the 
applications of these crystals. A number of studies 
have focused on experiment, but the constitutive 
model of ferroelectric single is absent. Huber[5] and 
Bhattacharya[6-7] et.al established a model based on 
the micromechanical method, but the simulation can 
not describe experimental result well. 

In this study, a polarization rotation model will be 
suggested to explain the stress-strain response of 
PMN-0.32PT. In order to validate that model, the 
stress-strain curve along <001> crystallographic 
direction of ferroelectric single crystals BaTi03 has 
been simulated in first principle method. Finally, 
based on the polarization rotation model, a 



constitutive model based on micromechanical method 
is constructed . 

Experimental Results and Discussion 

Experimental setup is shown in Fig. 1 to allow 
simultaneously imposing uniaxial stress and electric 
field to the specimen along the thickness direction. 
Mechanical load is applied by a servo-hydraulic materials 
test system (MTS). During test, the compressive stress 
loading and unloading cycles being -0.4— >-12— >- 
0.4MPa, -0.4->-40->-0.4MPa, and -0.4->-60->-0.4MPa, 
respectively. After each stress cycle, a unipolar electric 
field in the form of triangle wave with magnitude of 
0.5kV/mm and frequency of 0.02Hz is applied to the 
specimen to ensure the specimen return to its initial 
state. 
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FIG.l THE SKETCH of EXPERIMENTAL SET 
Measured a 33 - £ 33 responses of <001>-oriented PMN- 
0.32PT single crystal under the compressive stress 
loading and unloading cycles with the maximum 
magnitudes being 12MPa, 40MPa and 60MPa, 
respectively, are shown in Fig. 2. For the 40MPa and 
60MPa stress cycles (Fig. 2), the stress induced 
behavior of PMN-0.32PT single crystals shows shape 
memory feature. This can be explained by the 
polarization vector rotation mechanism. The detail 
description of polarization vector rotation model can 
refer to our paper [8]. 
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FIG.2 STRESS-STRAIN CURVES OF UNDER THE LOADING OF 12 
MPa,40MPa,60MPa 

First Principle Calculation of Stress-Strain 

Calculation Methodology 

PMN-PT and BaTi03 have the similar AB03 structure, 
so they have the similar ferroelectric properties. There 
is only Ti4+ particle in B site of BaTi03, however, there 
is not only Ti4+ but also minim Mn4+ and Ni2+ 
particle in B site of PMN-PT. So the single cell of 
BaTi03 is convenient in calculation and it keeps the 
similar ferroelectric to PMN-PT. 




[100] 



[010] 



FIG. 3 THE CELL CALCULAATION MODEL OF BaTi03 in R 
PHASE 

In this paper, we calculate a single cell of BaTi03, the 
single cell should be the smallest periodic reduplicate 
cell, it includes one Ti4+ particle, three 02- particles 
and one Ba2+ particle(Fig. 3). It is suggested that the 
initialized state of BaTi03 ferroelectric single crystal is 
R phase after {001} oriented polarization. Refer to 
literature[9], the crystal lattice constant of R phase 

BaTi03 ferroelectric single crystal is 4.001 A , the 
coordinate of particle in single cell is shown in Table 
l.The loading along {001} direction is carried out 
through application increasing strain by degrees. 
Stress and other parameters under each strain level is 
calculated by VASP. Calculation under each strain 
level include two steps. Firstly, the particle coordinate 
of single cell under each strain level is calculated by 
first principle molecular dynamic method. Secondly, 



Stress and other parameters are obtained through 
relaxation that is based on the first result. Each 
increment of strain in this paper is 0.5%, until 4%. 

TABLE 1 THE COORDINATE OF PARTICLE IN R PHASE BATI03 
FERROELECTRIC SINGLE CRYSTAL CELL 
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Discussion 

In order to validate that the method amentioned in 
this paper is correct, we calculate the elasticity 
constant C 33 of T phase BaTi03 using the method 
mentioned before. Fig. 4 shows that the C 33 is 180GPa, 
which is approach to the experimental value 
189 ± 8Gpa reported in refer[9]. This means that the 
method used in this paper is trusty. 





-0.5 
-1 

j? -1.5 

O -2 
$ -2.5 

-3 
-3.5 

-4 



-1 r- 



BaTiCX 



-2 -1.5 -1 -0.5 

£ 33 (%) 

FIG. 4 THE STRESS-STRAIN CURVE OF BATI03 FERROELECTRIC 
SINGLE CRYSTAL OBTAINED BY THE METHOD AMENTIONED 
IN THIS PAPER, THE {001} ORIENTED BATI03 FERROELECTRIC 
SINGLE CRYSTAL IS IN T PHASE 

Stress-strain curve of {001} orientated R phase BaTi03 
calculated in first principle method is shown in Fig. 
5(a). Fig. 5(a) shows that the result has similar nonlinear 
behavior to the experimental result of PMN-0.32PT 
during the loading that sketched in Fig. 2, namely, 
there is obvious "a,b,c" steps during loading. 
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FIG.5(A)The STRESS-STRAIN CURVE of (001} ORIENTED BaTi03 
FERROELECTRIC SINGLE CRYSTAL CALCULATED in FIRST 
PRINCIPLE CALCULATION, the SCATTER POINT is 
CALCULATION VALUE and the BROKEN LINE FIT the SCATTER 
POINT.(B) The SKETCH MAP of POLARIZATION ROTATION 
DURING the CALCULATION. The "A,B,C"LETTER in (B), the 
LETTER INDICATES DIFFERENT POLARIZATION STATES 
DURING LOADING. 

The Viscoplastic Model of Ferroelectirc Single 

The Viscoplastic Model 

From the experimental analysis, it is known that the 
<001>-oriented PMN-0.32PT single crystal undergo the 
R— M — T and R— M phase transform ation under 
the compression. In order to establish a compact 
model and keep the essence of experimerrtQ R 
phase transformation of PMN-0.32PT is considered in 
this study. Polarization rotation cause bulk 
deformation and shear deformation associated with 
the slip planes. It is assumed that the corresponding 
deformation associated with the slip planes is shear 
dominated, a feature similar to that of the multi-slip 
system of the crystal plasticity. This similarity renders 
possible to use crystal plasticity models to describe the 
transformation deformation of PMN-0.32PT single 
crystal. 

The poled ferroelectric single is four domain state, the 
polarization vector is along <111>,< 1 11>,<1 1 1> and 

< 1 1 1> respectively. In this study, it is assumed that 
the polarization vector of R phase can switch to the 
polarization vector of O phase, such as <111> vector 
switches to <001>, the vector also can switch back from 
O phase to R phase. The polarization rotation has an 
analogy to crystal plasticity slip, so we suggest there is 
eight slip systems in PMN-0.32PT ferroelectric single 
crystal. According to the crystalgraphic theory, we 
suggest ferroelectric single crystal has possible 8 
variants; transformation of variants can be 
characterized by the habit planes illustrated in Fig. 6. 
where n denotes the unit normal to the habit plane 
and s refers to the direction of transformation. 





FIG.6 THE COMPARISON BETWEEN FERROELECTRIC PHASE 
TRANSFORMATION AND PLASTIC SLIP OF SINGLE CRYSTAL. 
N IS DIRECTION NORMAL TO SLIP SURFACE. S IS DIRECTION 
ALONG PHASE TRANSFORMATION. 

In general, a criterion (i.e., the phase transformation 
criterion) exists for the phase transformation of 
ferroelectric single crystal, and the material is assumed 
to undergo phase transformation when at least one of 
the 8 variants satisfies the phase transformation 
criterion. This is detailed as follows. Upon loading, the 
condition to produce R^0 polarization rotation on a 
specified habit plane is that the driving force G of that 
plane reaches the critical value GOO . The driving force 
is composed of the chemical driving force Gchem and 
mechanical driving force Gmech 



G G chem + G mech G q 



(2) 



A similar condition holds for reverse transformation 
from CH-R polarization rotation with a critical value 
GOA 

G = G chem + G mech = G 0A ( 3 ) 

In Eqs. (2) and (3), the mechanical driving force can be 
expressed by 

Gi ch =Ty+E"P' (4) 

Where t" and E" are the resolved stress on the " a " 
transformation system, and y and P denote the 
associated transformation strain and transformation 
polarization. Following the crystal theory of plasticity, 

the resolved stress t" of the variant " a " is related to 
the stress tensor av. and the Schmid factor fl° by 



r a = V a a a 



(5) 



(6) 



i.7=l 

where the Schmid factor is defined as follow 

a ~\unn ■ mn ) 

with m i and n. being the unit normal to the habit 
plane and shear direction to the variant " a ", 

respectively. The resolved stress E" of the variant 
" a " is related to the electrical field tensor E, 



(7) 



The chemical driving force in Eqs. (2) and (3) is 
assumed to be a linear function of the temperature 
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where T and T denote the temperature in the single 
crystal and the equilibrium temperature respectively, 
and J3 is the stress-temperature coefficient. Note that 
the equilibrium temperature T Q is defined as the average 

of the starting temperature of O phase transformation 
and that of R phase transformation, that is 

r =i(o s +r) (9) 

When applying the rate independent crystal theory 
based model, Eqs. (2) and (3), to simulate the behavior 
of ferroelectric single crystals, one of the most 
computationally consuming tasks is to determine the 
set of instantaneously active transformation systems 
among the 8 possible variants at crystal level. This 
determination is usually achieved by an iterative 
procedure and must be carried out at each loading 
step, requiring extensive computation. Note that in the 
crystal theory of plasticity, a similar problem exists 
whilst in a rate dependent viscoplastic version of 
crystal theory of plasticity, determination of the set of 
active transformation systems is not necessary. As a 
result, computation effort can be reduced significantly. 
Following this idea, a viscoplastic version of Eqs. (10) 
and (11) are proposed and employed in this study. 
In the viscoplastic crystal model for PMN-0.32PT 
ferroelectric single crystals all transformation systems 
are assumed to be instantaneously active of varying 
extent, which is governed by a rate dependent viscoplastic 
law. In this paper, the phase transformation of variant 
" a " is assumed to comply with the following power 
law of viscoplasticity, 



r=/ 



G" 


— i 

" G" 


c 


G 


G 





(10) 



where / " is phase fraction transformation rate of variant 
" a ", f is reference phase fraction transformation rate, 
G" is the driving force of variant " a ", G is the 
critical driving force, exponents k and m are material 
parameters dictating the rate effect, c depends upon 
the phase fraction, c the reference value at initial 
state. 

The phase fraction of variant " a " is calculated as, 

r =/"//„ (11) 

Summation over all possible variants provides the 
phase fraction for the whole ferroelectric single crystal, 
that is 

£ = ±<r ( 12 ) 



The incremental form of Eq. (12) can be written as 

d^=df/f (13) 

The transformation strain tensor s". , which is associated 
with df , can be obtained as 

ds;=ta;fdr (14) 

r=l 

where a is Schmid factor, ds" is increment of phase 
train, y is the maximal strain transformation during 
loading, df" is increment of phase fraction. 

Elastic strain and electrical field induced strain during 
loading can expressed as 



ds e . =C..,,da ... ds -d, dE, 

ij ijkl ij ' ij kij k 



(15) 



Increment of the total strain ds v can then expressed as 
the sum of an elastic component ds'.. , a electrical field 
induced component ds^ and a transformation induced 
component ds" , 

ds = ds* + ds E + ds" (16) 

ij ij ij ij \ / 

Numerical Results 

To validate the model, the corresponding calculation 
result is compared to experimental stress-strain curve. 
The material parameters used in this study are 
as follow[10,ll]. Material parameters of R phase 
are C* =9.2GPa , C^ 2 =10.3GPa , C 4 * 4 =6.9GPa , 
y = 0.0033 , G 0R = 0.5073MJ/m 3 , /? = 0.004MPa/K , 
l/m = ll , l/jt = 1.0 . T = 403 , J =298K . Material 
parameters of O phase are C° = 38GPa , C° = 40GPa , 
C 4 ° 4 = 28GPa , G 0o = 0.276MJ/m 3 , = 0.039MPa/K , 

l/m = 11.5 , l/k=l.l . Fig. 12 shows the model 
predicted and experimental measured response. 
Result in Fig. 7 shows that stress-strain response of 
PMN-0.32PT can be predicted by the developed 
constitutive model, with quantitative agreement. 
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FIG.7 THE COMPARISON BETWEEN EXPERIMENTAL AND 
SIMULATION STRESS-STRAIN CURVE ALONG <001> 
DIRECTION OF FERROELECTRIC SINGLE CRYSTAL. 
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Conclusions 

In this paper, polarization rotation model is developed 
to explain the observed behaviors of PMN-0.32PT. The 
stress-strain curves along <001> crystallographic 
direction of ferroelectric single crystal BaTiOs is 
calculated with the first principle method. Obtained 
results show that the R— »M— »0 polarization rotation 
takes place in BaTiCb ferroelectric single crystals under 
compression, which is consistent with the polarization 
rotation model. Based on the polarization rotation 
model, a constitutive model based on micro-mechanical 
model is proposed. It is shown that the developed 
model can faithfully capture the typical characteristic 
of the observed constitutive behavior of <001>- 
oriented PMN-0.32PT 
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